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Abstract. In this paper we calculate a number of integrals over SU(iV) of interest 

■ in Hamiltonian Lattice Gauge Theory calculations. 
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1. Introduction: file preparation and submission 

This paper is concerned with the analytic calculation of certain group integrals over 
SU(iV) which are of interest in Hamiltonian Lattice Gauge Theory (LGT) pQ. In par- 
ticular, the integrals in question arise in calculations of ground state energies and mass 
gaps using the plaquette expansion method [2] . While these integrals have traditionally 
been handled with with numerical Monte Carlo integration, in this paper we demon- 
strate how these integrals may be calculated analytically. 

Traditionally Hamiltonian LGT techniques have been restricted to simple "trial" 
states (most commonly the strong coupling ground state which is anihilated by the 
chromo-electric field). The techniques presented in this paper are are of use in the ex- 
tension of standard Hamiltonian LGT techniques beyond trivial "trial" states, at least 
in 2+1 dimensions, to include more complicated trial states than the stong coupling 
ground state. For example, the techniques presented in Section [2] have been applied 
in variational estimates of glueball masses in 2+1 dimensions [31 H]. It is thought that 
similar techniques may be of use in plaquette expansion method calculations of glueball 
masses in 2+1 dimensions. 

Much work has been carried out on the topic of integration over the classical com- 
pact groups. The subject has been studied in great depth in the context of random 
matrices and combinatorics. Many analytic results in terms of determinants are avail- 
able for integrals of various functions over unitary, orthogonal and symplectic groups [5]. 
Unfortunately similar results for SU(iV) are not to our knowledge available. Their pri- 
mary use has been in the study of Ulam's problem concerning the distribution of the 
length of the longest increasing subsequence in permutation groups [6j [7] . Connections 
between random permutations and Young tableaux [8] allow an interesting approach 
to combinatorial problems involving Young tableaux. A problem of particular inter- 
est is the counting of Young tableaux of bounded height [9] which is closely related to 
the problem of counting singlets in product representations. Group integrals similar to 
those examined in this paper have also appeared in studies of the distributions of the 
eigenvalues of random matrices [TO], [TT] . 

In the last 20 years it appears that not much effort has been devoted to the subject 
of group integration in the context of LGT. The last significant development was due 
to Creutz who developed a diagrammatic technique for calculating specific SU(iV) inte- 
grals p2] using link variables. This technique allows strong coupling matrix elements to 
be calculated for SU(iV) [13] and has more recently been used in the loop formulation 
of quantum gravity where spin networks are of interest [HI [151 03] • 

This paper is structured as follows. Section [2] contains calculations of some general 
SU(iV) integrals which are used in later sections. Section [3] contains integrals of SU(iV) 
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group elements which are then applied in Section H] to the calculation of a number of 
integrals that arise in the plaquette expansion method in 2+1 dimensional Hamiltonian 
LGT using a one plaquette trial state. 

2. Preliminaries 

In an earlier paper [3] a useful technique for performing SU(iV) integrals was described. 
In this paper the following results in terms of Toeplitz determinants were derived, 

Gxj (N) (c,d) = I dUe cTrU+dTrU ' = det [/,-_< f 2Vcd)l (1) 

JV(N) l J \ /J l<i,j<N 

Gsu(iv) (c, d) = / dUe<™^ = V - det ( 2 VS) 1 . (2) 

Here c is any number, dll is the Haar measure normalised so that J dU = 1, the 
quantities inside the determinant are to be interpreted as the (i, j)-th entry of an iV x iV 
matrix, and /„ is the n-th order modified Bessel function of the first kind defined, for 
integers n, by 

00 x 2k+n 

Eq. (j2J) can be used to as a generating function to calculate integrals of the form 
/ (TiU) n {TiU ] ) m dUe cTtU+dTrU \ (4) 

JSV(N) 

by direct differentiation. This process allows easy calculation of integrals of polynomials 
in TrU and TrW. It does not allow the simple calculation of integrals involving the more 
complicated trace variables, (Tr?7 n ) m for integers n and m. 

Making use of the same technique described in J3], the following more generals 
integrals can be calculated, 

dU Xni n 2 -n N {U)e cTTU+d ™' = ( - J det Ij_ i+ni (2 V^d) and (5) 

U(A0 



c / 

/ dU Xni „(Uy™ +d ™' = - Yl - det/ i _ w+rii (2V^). (6) 

To arrive at these results, instead of Eq. (15) in [I], we make use of the following result 
due to Weyl [17] (with implicit sums over repeated indices understood) and follow the 
same procedure, 

A(0i, . . . , (j> N )Xn,-n N = - ? L=£ ilJ2 ... ijV e^( iV - il+ ^)e^( iV -^ + ^) ■ ■ ■ e ^(N-i N+ n lN )_ (?) 

Here e^...^ is the totally antisymmetric Levi-Civita tensor defined to be 1 if {ii, ... , i n } 
is an even permutation of {1,2, ... ,n}, —1 if it is an odd permutation and otherwise 
(i.e. if an index is repeated). 
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Eqs. and allow the simple calculation of many integrals involving (TiU n ) m 
for integers n and m. To perform these calculations expressions of the form (Tr£/ n ) m 
must be written in terms of characters. To do this the following results are useful 

mm, 

Xnr a ...r lf - 1 {U) = det [xri-H-i]i<ij<jv-i , where (8) 
Xn(U) = J: 5 I Yikt - n ) f[ j^iWp. (9) 



/ n \ n 1 



Using these results, expressions for characters in terms of trace variables can be 
calculated. These results can be rearranged to give expressions for trace variables in 
terms of characters. As an example, for SU(iV) it can be shown that 

(Trf/) 2 = X2{U) + X u(U) and (10) 
Tt(U 2 )= X 2(U)-Xii(U). (11) 
Making use of Eqs. (|T0|) . (TTTT) and the following integrals can be easily calculated 

/ dUTr(U 2 )e cTrU+dTvUt and / dU(TrU) 2 e criiU+dTlU \ (12) 
3. Creutz's Method 

In this section we make use of Creutz's method for calculating strong coupling 
integrals over SU(iV) to calculate more complicated integrals. In the terminology 
of Hamiltonian LGT, we use Creutz's method to calculate one-plaquette trial state 
integrals. The integrals we consider in this section are of use in the calculation of matrix 
elements arising in Hamiltonian LGT, particularly in calculations involving Hamiltonian 
moments. 

The first two integrals we consider calculate in this section are, 

W*(c) = J dUe c ^ u+u ^U ab U a , v and (13) 

J ab -Mc) = I dUe cT *V +u ^U ab (rf) a , b ,, (14) 



where c is a complex valued variable. 

Let us start with I a b;a'b'(c)- Expanding the exponential in a Taylor series and also 
expanding the resulting binomial terms gives, 



oo m I \ m f 

Wf(c) = E E 7 )~\ / dUUaiai ' ' ' U *~*~iW)*& ■ ■ ■ ( U \kUa b U a > b >. (15) 

m=0 i=0 \ / ' ^ 

Following the diagrammatic procedure first described by Creutz [TJ], the contributing 
integrals within the sum reduce to products of paired Levi-Civita symbols. The indices 
of each pair of Levi-Civita symbols are completely contracted with one another with 
the exception of those that contain an uncontracted U index a, b, a' or b'. Following 
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Creutz's procedure, each non-zero integral under the sum in Eq. (TT5"j) can be reduced to 
linear combinations of the following products of Levi-Civita symbols, 

£aA 1 A 2 ---A N _ 1 £bA 1 A2---A N _ 1 £a>B 1 B2---B N _ 1 £b'B 1 B 2 ---B N - 1 ^ 4ib<5a'fc' , arid (16) 
£aa'A 1 A2---A N _2 e bb'A 1 A 2 ---A N ~2 ^ &abO~a'b' ~ 0~ab'0~ba', (17) 

as the paired Levi-Civita symbols can contract with uncontracted U index pairs, (a, b) 
and (a', b') in two possible ways, one involving a single pair of Levi-Civita symbols 
(Eq. ( 1161) ) and the other involving two (Eq. (1T7I) ). The additional, fully contracted 
Levi-Civita symbols contrbute a multiplicative factor to each term. The simplification 
to delta functions presented in Eqs. (TlBl and ( fTTl) is easily obtained from the standard 
result from differential geometry, 

Hence, I a b;a'b> may be expressed as a linear combination of delta functions, 

Iab-,a'b'(c) = bi(c)5 a b5 a 'b' + b 2 (c)5 ab > 5 ba > , (19) 

where &i and 62 are functions of c to be determined. To find these functions we contract 
over indices in Eq. ( fl9l) to contruct equations for b\ and b 2 in terms of known integrals. 
For this particular integral we construct two equations for b± and b 2 by firstly, mulitpling 
both sides of Eq. ( fl9l) by 8 ba 8y a i and summing over repeated indices and secondly, 
mulitplying both sides of Eq. ( !T9l by S a ' b S b ' a and summing over repeated indices. These 
operations result in the following two equations: 

Iaa;bb(c) = N%( C ) + Nb 2 ( C ) (20) 
Iab;ba(c) = Nb 1 (c)+N 2 b 2 ( C ) J (21) 

which can be solved for b\ and b 2 and substituted in Eq. ( fl9l to give our final result for 
Iab;a'b' m terms of known integrals, 

Tr(U 2 



L ab:a'b' 



( C ) = | d[/e cTt ( c/+c/t ) 

+ AT2_1 J aUe 



(TrUy - 
Tr(U 2 



N 
(TrU) 



(22) 
(23) 



N 

Our approach to calculating J a b;a'b', defined in Eq. flHj) , follows the approach for 
I a b;a'b> decribed above. As was done for I a b;a'b', it can be shown that Jab-a'V m ay be 
expressed as a linear combination of delta functions using the method due to Creutz, 

Jab;a'b'{c) = di(c) 6 ab 5 ' a > y + d 2 (c)5 abl 5 ba i . (24) 

By contracting with appropriate delta functions we arrive at simultaneous equations for 
d\ and d 2 , which are solved to give an expression for J ab] a'b' in terms of known integrals 

Jab-Mc) = / dUe c ^ u+u ^ (TrC/TrC/t - l) (25) 

+ [dUe cT ^ u+u ^ (n- ±-TrUTrtf) . (26) 



N 2 - 
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This result agrees with the strong coupling result J a fe ; a'6'(0) = o~ a b>5b a i/N since 



dUTrUTrU* 



dU = 1. 



(27) 



This procedure can be used to evaluate more complicated integrals. For example, 
following Creutz's method again, 



Iab;a'b';a"b"(c) — J dUe cTl< " U+U ^U a bU a 'b'U a "b" , 



(28) 



can be written as as linear combinations of the following products of Levi-Civita symbols, 

£aA 1 - A N - 1 £bA 1 -A N _ 1 £a'B 1 - B N - 1 £b>B 1 - B N - 1 £a"C 1 -C N - 1 £b"C 1 -C N _ 1 OC S a bS a 'b'5 a »b" , (29) 

£aa'A 1 ---A N -2 € bb'A 1 ---A N ^2 e a."B 1 ---B N ^ 1 £b"B 1 ---B N _ 1 OC {b~ a bO~a'b' ~ 0~ab> ' &ba')o~a"b" , (30) 

£aa"A 1 -A N „2 e bb"A 1 -A N „2 e a'B 1 -B N ^ 1 £b'B 1 -B N ^ 1 OC (b~ a bb~a"b" — 5 a b"5 b a")o~a'b' , (31) 

£a'a"A 1 -A N -2 e b'b"A 1 -A N „2 e aB 1 -B N _ 1 £bB 1 -B N „ 1 OC (S a 'b'S a " b 'i — S a > b "5 b ' a ")5 ab , and (32) 
^aa'a"Ai-A N ^bb'b"Ai-A N - 3 OC 8 ab " (5 a 'b>8 a " b ~ 0~a'bO~a"b') + 0~ab' (~ 8 a 'b"0~a"b + 0~a'bO~a"b") 

+8ab(8a'b"0~a"b' — 0~a'b'0~a"b")- (33) 

These results are easily derived using Eq. (fl8"|) . Hence, I a b;a'b';a"b" m ay be written as the 
following linear combination 

Iab;a'b';a"b"(c) = X\8 a b" 8 a' V $ a" b + X 2 8 ab "8 a 'b8 a " b i + X36 a b'8a'b"8a,"b 

+ X4d~ a bi5 a 'b0~ a "b" + X^8 a b8a'b"8a"b' + XQ8 a b8a'b'8a,"b" ■ (34) 

A collection of six equations for the Xi coefficients may be constructed by contracting 
appropriate indices as follows: 

/ N 2 N N N 2 N 2 N 3 \ 

N N 2 N 2 N N 3 N 2 

N 3 N 2 N 2 N N N 2 

N N 2 N 2 N 3 N N 2 

N 2 N 3 N N 2 N 2 N 

\ N 2 N N 3 N 2 N 2 N ) 

This system may be solved to give, 

2i 2 + 2i 4 - N(h +i 5 + i 6 ) + (N 2 - 2)z 3 



X-2 

X4 
X 5 

\x 6 J 



^ Iaa;bb;cc ^ 
Iaa;bc;cb 
Iab;cc;ba 
Iab;ba;cc 
Iab;bc;ca 
\ Iab\ca\bc ) 



Ik \ 

k 
k 
k 



(35) 



x x 
x 2 
x 3 

x 5 
x 6 





(N- 


2)(N- 


1)N(N + l)(JV + 2) 




2ii 


+ 2i 6 


-N(i 2 


+ iz + i 4 ) + {N 2 -2\ 


k 




(N- 


2)(N- 


1)N(N + 1)(N + 2) 




2zi 


+ 2i 5 


-N(t 2 + t 3 + t 4 ) + (N 2 -2] 


k 




(N- 


2)(N- 


l)JV(iV + l)(JV + 2) 




2i 2 


+ 2i 3 


-n(h 


+ l5 + i6 ) + (N 2 -2] 


k 




(N- 


2)(N- 


1)N(N + l)(JV + 2) 




2i 3 


+ 2i 4 


-N(h 


+ ^ + ^ 6 ) + (iv 2 -2; 


k 




(N- 


2)(N- 


1)N(N + 1)(N + 2) 




2i 5 


+ 2i 6 


-N(i 2 + t 3 + t 4 ) + (N 2 -2] 


k 




(N- 


2)(N- 


1)N(N + l)(JV + 2) 





(36) 
(37) 
(38) 
(39) 
(40) 
(41) 
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Making use of the definitions in Eq. (1351) . we see that i 2 = i 3 = 14 and i 5 = i e . Making 
these replacements in Eqs. fl36l— [4"T]) and substituting the results in Eq. fl34|) gives, 
T (A -iV(^i + 2z 5 ) + (iV 2 + 2)z 2 ^ 

1 ab;a'b';a"b"[ c ) ~ _/V(_/V"2 1)(]V 2 4) + ab 1 a' b<J a" b" + a b a* b" a" b> ) 

2i 1 -3Ni 2 + N\ 

+ N(N 2 - 1)(N 2 - 4) l daft " da ' ftda " 6 ' + d a&'<W0o"6j 

4^ 5 - 3iV» 2 + (iV 2 - 2)1! 
+ iV(iV 2 -l)(iV 2 -4) < U < W < W - ( 42 > 
The integrals, i 2 and i$ can all be calculated using the procedure set out in Section [2j 

In order to calculate all integrals of three Uij or (U^)ki factors only one more integral 
is required, 

J ab ;a>V;a»b»(c) = [ dUe C ^ U+U ^U ab (U^) a r V U a , ¥ r . (43) 



Other integrals of three group elements can be calculated using results for I a b;a'b';a"b" 
and J a b;a'b';a"b" by making appropriate changes of variables. 

J ab ;a'b';a"b" can be calculated in the same way as I a b;a>b';a" b ", in fact J a b;a'b>;a"b" is a 
linear combination of the same delta functions as I ab ;a'b';a"b" with coefficients defined 
by Eqs. (1361 T4T]) but with the definitions of i%, . . . , z 6 modified to reflect the integral in 
question (i.e. i\ = J aa ;bb,cc, ■ ■ ■ = J a b,ca;ba) ■ The final result can be shown to be 
2Ni 5 - Ni x + (N 2 - 2)i 3 
N(N 2 - 1){N 2 -4) 
2ix - Ni 3 - N\ 



Jab;a'b';a"b"(c) — AT/ Aro ; N , ; TO — 0~ ab" 0~ a' V $ a" 



N(N 2 - 1)(N 2 -4) 
2z 3 - JVii + N(N 2 - 2)25 



(<W'5 a '6(5a"&' + 0~ab'0~a'b"0~a"b) 

(0~ab'0~a'b0~a"b" + & ab& a' b" & a" b') 



N(N 2 - 1)(N 2 -4) 
2( 2 - N 2 )i, - jVz 3 + (iV 2 - 2)h 
+ iV(iV 2 -l)(iV 2 -4) d a*<WW. ( 44 ) 

Using a software package cabable of symbolic manipulation, such as Mathematica, this 
technique can be easily extended to more complicated integrals, involving four or more 
U and / or TP group elements in the integrand. In order to accelerate the calculation 
of such integrals, /, the following procedure is used: 

(i) Express linear combination of appropriate delta function products, 
with unknown coefficients, Xj. 

(ii) Simplify this expression, taking into account symmetries of the integral with respect 
to {ax, • • • , a m } and {bx, ■ ■ ■ , b m }. Some Xi will be found to be equal to others, 
reducing the number of independent Xi. For example, from Eq. ( 1281) . we see that 
Iab;a'b';a"b" = 1 a" b" ;a' v ;ab ■ Making use of this and Eq. (1341) we find that x 2 = x 3 and 
Xx = x 5 . 

(iii) Construct a system of equations for x^ in terms of known integrals by contracting 
appropriate combinations of indices. 
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(iv) Solve this system for Xi and substitute into the original expression for I i6i;--;a m 6 m 
from step @. 



4. Integrals Required for Moment Calculations 

In this section we calculate a number of integrals useful in Hamiltonian LGT, especially 
for those techniques requiring the calculation of Hamiltonian moments. The calculations 
make use of the results of Section [31 



The Hamiltonian methods for LGT generally require the calculation of various in- 
tegrals over SU(iV). The t-expansion [19] and plaquette expansion [2] methods are two 
approaches to Hamiltonian LGT which require the calculation of Hamiltonian moments. 
These moments can be simplified so that all electric field operators are removed, and 
thus reduced to a set of group integrals over SU(iV). For simple one-plaquette trial states 
in 2+1 dimensions these integrals can be calculated analytically. More complicated sce- 
narios (3+1 dimensions or more complicated trial stated) typically require numerical 
evaluation. 



In the calculation of Hamiltonian moments in 2+1 dimensions, the following 
integrals require calculation at the fourth order for diagrams on two-plaquette skeletons 
(in the following U and V represent the neighbouring plaquettes of the two-plaquette 
connected skeleton): 



Z 1 (c,N) =<0d 

Z 2 (c,N) =(</>o 
Z 3 (c,N) 

Z 4 (c,N) =(0 O 

Z 5 (c,N) =(0 O 

Z 6 (c,N) =(0 O 

Z 7 (c,N) =(0 O 

Z 8 {c,N) =(<P' 

Z 9 {c,N) =<0' o 



Tr(UV)Tr(UV)\(j)' ) 



(0o|0o) 



Tc(UV)Tt(UV^M) = ( -DO- > (0 o |0o) 



Ti{UV ] )Ti{UV ] )\(f)' Q ) 
Tr^^Tr^V^l^d) 



(0o|0o) 



(0o|0o) 



TriUUVVM) e 
Tr(UUV j V j )\(f)' } = 



Tr(UVUVM) = 



Tr(UV^UV)\(f>' Q ) = 
Tr(UVWV)\((>' ) = 



(0o 1 00 ) 



(0o|0o) 



(00 1 00 ) 



(00 1 00 ) 



Wo/ 



(45) 
(46) 
(47) 
(48) 
(49) 
(50) 
(51) 
(52) 
(53) 
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Z 10 (c,N) = (<f>' \Tr{UV^UV^M) 



C|3} ) (<h\<M 



(54) 



Here each integral is calculated over SU(iV) and |0o) is the one plaquette trial state: 



)) = exp 



10), 



(55) 



where |0) is the strong coupling vacuum and the sum is over all plaquettes on the lattice. 
|0q) is the reduced trial state incorporating only the neighbouring plaquettes U and V, 



!<#,)= exp ^-Tt(U + E/+ + V + V^) |0). 



L2 



(56) 



Eqs. (SSJ), f|4T|) . (jUD , (J5l~j) and fl54l) can be calculated using standard character 
expansion techniques pE]. Here we will calculate these integrals using a different method 
which is able to be applied to the calculation of the other integrals presented above. 
We start with Zi(c,N), which can expressed as (assuming summation over repeated 
indices) : 

= Iab;a'b'(c)Iba;b'a'(c). (57) 

Making use of Eqs. (TT9]) and (1231) . Z\ can be reduced to an expression involving known 
single plaquette integrals as follows: 

Zl(c,N) = (Ma6<W + &2<W<W)(Mba<W + b 2 5 ba '5 a y) 

= N 2 b\ + 2N 0l b 2 + N 2 b\ 



NILm 



2I aa ;bb^ab;ba :ba 



(58) 
(59) 

(60) 



Z 2 (c,N) 



(61) 
(62) 



N(N 2 - 1) 

Next we calculate Z 2 . Using the same technique we have, 

' dUdVe^ u+u \^ v+v ^U ab V ba UMV%^ 

= Iab;a>b'(c)Jba;b'a'(c). 

Making use of Eqs. (fT9|) and (1231) and also Eqs. (124")) and (!26l) . Z 2 can be reduced to an 
expression involving known integrals as follows, 

Z 2 (C, N) = Ioh-a'V (c) Jba;b'a' (c) (63) 

= iV 2 Mi + + iVMi + N%d 2 (64) 

N I aa\bbJ aa\bb I aa\bbJ ab\ba Iab:ba,Jaa;bb N 1 ab\baJ ab\ba 



N(N 2 - 1) 



(65) 



Z 3 can be shown to be equal to Z\ as follows, 



Idc 



Z 3 (c,N) = J dUdVe^ v ^e^ v ^U^{V% a {V\ 

= Iab;cd(c)Iba;dc(c) = Zx(c,N). (66) 

Here the last line follows from the fact that an integral of a single plaquette loop does 
not depend on its direction. 
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The remainder of the integrals can be calculated similarly. The results are as 
follows: 

1 2 

Z 4 (c, AT) = J^——^(J aa . bb + J ab . ba ) — _ Jab;ba,Jaa;bb (67) 

Z 5 (c,iV) =^4, a (68) 

Z 6 (c,N) =Z 5 (c,N) (69) 

9 NT T _ — 

^ / ArN z,iv 1 aa;bb-l ab;ba 1 aa;bb 1 ab:ba 

Z ^ N) = N{N^T) (70) 

ry i .r\ N I aa.bb Jab.ba I aa;bbJ aa;bb I abibaJ ab\ba N I abba J aabb t<-ii\ 

Z§(C ' N) = N(N^T) (71) 

Z 9 (c, N) = 2iYJaa;febJ ^;5;- feb " JlbM (72) 

Z 10 {c,N) = Z 6 {c,N) (73) 

The techniques presented in this section can be extended to the calculation of loop 

integrals on three or more neighbouring plaquettes using Eqs. ( H2|) and (jUJ) and their 
extensions. 



5. Conclusion 



In this paper a number of integrals over SU(iV) have been calculated. Further work 
could incorporate these calculations into plaquette expansion calculations of glueball 
masses in 2+1 dimensions. Other calculations involving Hamiltonian moments may be 
able to be extended to more complicated trial states, at least in 2+1 dimensions, using 
the techniques presented in this paper. 
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